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Abstract. A thermodynamic analogue of the Pauli problem (re- 
construction of a wavefunction from the position and momentum 
distributions) is formulated. The coordinates of a quantum sys- 
tem are replaced by the inverse absolute temperature and other 
intensive quantities, and the Planck constant is replaced by the 
Boltzmann constant multiplied by two. A new natural mathemat- 
ical generalization of the quasithermodynamic fluctuation theory is 
suggested and sufficient conditions for the existence of asymptotic 
<— j . solutions of the thermodynamic Pauli problem are obtained. 



1. Introduction 

Anyone who had studied statistical thermodynamics could have no- 
ticed a certain analogy between the mathematical formulae describing 
the quasithermodynamic fluctuations and the quantum mechanical for- 
mulae related related to coherent states and the Heisenberg uncertainty 
principle. Take a simple example of an abstract thermodynamic system 
with an entropy function 

S = S(U,V,N), 

where U is internal energy, V is volume, and N is the number of moles 
of a chemical substance (one may think of a Van der Waals gas). If 
we fix the volume V = Vq and the number of moles N = N Q , and put 
the system in a thermostat with absolute temperature T = T , then 
the internal energy U is going to experience fluctuations 8U around 
a value U = Uq satisfying an equation Tq" 1 = dS(U,Vo, N )/dU (we 
assume this solution is unique). The corresponding probability density 
fsu(x) in a point x G K. is approximately described by the normal 
distribution (Einstein's formula): 

fsu(x) = (2irk B [~S^ u ]- 1 )- 1 / 2 exp(S^ u x 2 /(2k B )), 

where k B is the Boltzmann constant, and S'^jj is a short notation for 
the second derivative d 2 S(U, V, N)/dU 2 taken in U = U , V = V , 
N = N (we assume that S'^u < 0). The variance of 5U is then as 
follows: 

{{5U) 2 ) = k B [-S'>v\-\ 
1 
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It is quite important to stress that the internal energy U is not the 
only quantity that is going to experience fluctuations. We have at 
the same time a fluctuation 5(3 of the inverse absolute temperature 
(3 = dS(U,V , N )/dU around the value (3 := l/T . Its probability 
density fsp{y) in a point y G K is approximated by the formula: 

My) = (2nk B [-S^ u ])- 1 / 2 exp(y 2 /(2k B S'^ u )). 

Note that fsp(y) coincides with the probability density of the linearly 
transformed random variable S'^jjSU. For the variance of 5(3 we have: 

{(5(3f) = kBl-Sfa], 

so we obtain: 

vf^vff) = k B . 
All this is very similar to the Heisenberg uncertainty relation in a 
coherent state. This observation was made by many authors, see for 
example [H [2], [3], SI El El . In fact the discussion of this analogy 
can be traced back to W. Heisenberg and N.Bohr (the complementarity 
principle). 

Take an abstract one-dimensional quantum mechanical system de- 
scribed by a wavefunction 




where a > is a parameter, q varies over R corresponding to the possi- 
ble values of a coordinate Q, ft is the Planck constant, and i\>% G L 2 (R). 
The self-adjoint operator Q representing Q acts as a multiplication by q, 
and the self-adjoint operator P representing the canonically conjugate 
momentum P acts as a derivation —ihd/dq. The state corresponding 
to ipn is coherent, i.e. we reach the lower bound in the Heisenberg 
uncertainty relation: 

^WQWnVWPWn = ft/2, 

where (—}n denotes the quantum mechanical average, 5Q = Q — (Q)h, 
and 5P = P- (P) h . 

The analogy between the two formulae is too striking to neglect. The 
Planck constant ft corresponds to the Boltzmann constant k B multi- 
plied by two, but to develop this analogy one needs answer the fol- 
lowing question: What is a "thermodynamic wavefunction"'! In the 
present paper I suggest to attack this question via the so-called Pauli 
problem. The Pauli problem is basically a problem of a reconstruction 
of a quantum mechanical wave function based on the directly available 
experimental data. Around 1933 W. Pauli has asked the following: 
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to what extent is a quantum mechanical wavefunction ip(q) G L 2 (R) 
determined by the corresponding marginal distributions 

I(q) = mq)\ 2 , j(p) = m P )\ 2 , 

of coordinates q G R and momenta p G R, 

The probability densities I(q) and J(p) (assuming ||^|| = 1) is some- 
thing we can see directly in an experiment (modulo the subsets of 
Lebesgue measure zero). Therefore the task is to reconstruct a wave- 
function ip from a given pair of probability densities {I, J). It turns 
out that the problem is quite complicated (W. Pauli himself had not 
given a complete answer). A generalization for a <i-dimensional quan- 
tum system is straightforward. It is also of interest to point out that 
there exists a tomographic generalization of this problem (for a review, 
see j8]): one is required to reconstruct a density operator p on L 2 {R), 
pi" = p ^ 0, Trp = 1, rather than a wavefunction ip, based on the 
marginal densities I^ u (z), z G R, corresponding to the observables 
p,Q + uP, where p and v vary over R. This problem is known to have 
a solution. 

It turns out that a solution of the Pauli problem to construct tp out 
of a given pair (J, J) does not need to be unique. Take as an example 
the following function: 

i>a,b,c(q) = (2a/7r) 1/4 exp(-2(a + ib)q 2 + icq), 
where a G R>o, b G R\{0}, c G R. It is straightforward to check that 

|t/V>,c(<?)| 2 = |^a,-6, C (g)| 2 , \^a,bA P )\ 2 = IVV-6,c(p)| 2 , 

for q,p G R. 

The main idea of the present paper is to consider the experimental 
data about the fluctuations of internal energy 5U and inverse absolute 
temperature 5(3 as an input data for the Pauli problem. More precisely, 
we may assume that we are given a pair of probability densities wsu{x) 
and w$p{y) of a more general shape than the Gauss exponents fsu{ x ) 
and fsp(y). The aim is to span a "wavefunction" over wsu(%) an d 
wsp(y). In the main part of the paper I investigate when this is actually 
possible, i.e. when does a thermodynamic wavefunction exist. It turns 
out that typically such a function can not be constructed, but on the 
other hand this only implies that the corresponding state is not pure. 
In other words, for a complete picture one needs to investigate a d- 
dimensional quantum tomographic problem in thermodynamics. 
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On the other hand, a thermodynamic wavefunction exists if we as- 
sume that the fluctuations are described by the Gauss exponents fsu{ x ) 
and fsp(y) precisely. Furthermore, it is known that the quasithermo- 
dynamic fluctuation theory based on these exponents is quite useful 
in practice and really works. Therefore it is natural to consider an 
analogue of a semiclassical asymptotics H — > which in this case we 
may denote as 2ks —> 0. In the second part of the paper I study the 
asymptotics of the thermodynamic Pauli problem corresponding to a 
vanishing Boltzmann constant and formulate a sufficient condition of 
existence of asymptotic solutions. The non-uniqueness of a solution in 
this case yields, in particular, a new mathematical generalization of the 
quasithermodynamic fluctuation theory. 

The general motivation for the present work comes from the idea 
to "quantize" thermodynamics in the sense of V.P.Maslov [91 [101 EH 
[121 [13]. Note that in [H] the "thermodynamic" terminology (enthalpy, 
entropy, etc.) is used to describe a (non-linear) complex germ on a 
Lagrangian manifold. One can look at the Lagrangian manifolds which 
emerge in thermodynamics and apply the semiclassical methods (the 
tunnel canonical operator [T3]) to study the asymptotics of the partition 
function of a system with respect to a large number of particles. A 
generalization of thermodynamics in this spirit is also of great interest 
in the context of quantum gravity [TH [T71 [THJ [19] . 

Note that "quantization" of thermodynamics uses a different quan- 
tization parameter than H. In this sense it is not exactly "quantum" 
thermodynamics (the thermodynamics of small systems) where the aim 
is to construct a kind of mixture (a grand unification) of quantum me- 
chanical and thermodynamic pictures J2U1 ED E21 E21 EI] (i-e. both h 
and are involved). At the same time, a thermodynamic wavefunc- 
tion in the sense as considered in the present paper (a "deformation" 
of quasithermodynamics) is an extremely natural concept. 

2. Thermodynamic Lagrangian manifolds 

The concept of a Lagrangian manifold [21 [TJ] provides perhaps the 
most natural way to axiomatize phenomenological thermodynamics. 
Note that in phenomenological thermodynamics we do not have a no- 
tion of a "number of particles" and the Boltzmann constant ks is also 
not there. Instead, we work with the number of moles N of a chemi- 
cal substance. In this sense, phenomenological thermodynamics is re- 
lated to statistical thermodynamics like classical mechanics is related to 
quantum mechanics. The importance of this analogy is stressed in [9]. 
The idea of "quantization" of energy (i.e. partitioning it into quanta) 
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leads to an introduction of the Planck constant h. Similarly, the idea of 
"quantization" of matter (i.e. partitioning it into corpuscles) leads to 
an introduction of the Boltzmann constant k B - The number of particles 
(corpuscles) in one mole is the Avogadro number, 

N A = R/k B , 

where R is the universal gas constant - a physical constant present at 
the level of phenomenological thermodynamics. 

Take a simple one-component thermodynamic system described by 
absolute temperature T, volume V, and number of moles N. Denote: 
U - the internal energy, S - the entropy, p - the pressure, and p - the 
chemical potential. Then on the states of thermodynamic equilibrium 
we have: 

dU = TdS -pdV + pdN. 
It is convenient to denote 

1_ ^__p_ „ _ -p 

P ■— rpl P -~ rpl t 1 - _ rp ) 

and to rewrite this equation as 

dS = f3dU + pdV + pdN. 

Assume that (U, V, N) varies over a domain D C IR 3 (£/, V, N) and that 
S is described as a smooth function S = S(U,V, N). Then we ob- 
tain a 3-dimensional surface A C R 6 (/3,p, p,, U, V, N) described by the 
equations: 

dS(U, V, N) dS(U, V,N) ^ _ dS(U, V, N) 

1 du ' p ~ dv ' /i ~ m • 

It is a common practice to write these derivatives as (dSfdU)v.N, 
{dS/dV)u.N, an d (dS/dN)u. v , respectively. The surface A is a La- 
grangian manifold with respect to the symplectic structure 

u := df3 AdU + dpAdV + dpAdN, 

i.e. the pull-back i* A (u) = 0, where i A : A ->■ R 6 (/3,p, p, U, V, N) is the 
canonical embedding. The first law of thermodynamics says basically 
that the differential 1-form a on A, 

a :=i* A (/3dU + pdV + pdN), 

is exact. The entropy S is an action on this Lagrangian manifold A 
and it is measured in the same units as ks- 



[S] = [k B ]. 
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Note that this is totally similar to mechanics: the mechanical action A 
of a system is measured in the same units as h, 

[A] = [h], 

and if we write A as a function of coordinates q — (qo, qi, ■ ■ ■ , q n -i) £ R™ 
(where n is the number of degrees of freedom), A = A(q), then we 
obtain a Lagrangian manifold L C M? n (p,q), p = (p ,Pi, • • • ,Pn-i), 
described by the equations pi = dA(q)/dqi, i = 0, 1, ... ,n — 1 (the 
symplectic structure is given by Y17=o d Pi A daft 

Let us define an abstract thermodynamic system as follows. Consider 
a phase space ]R 2 ( d+1 )(ft E), f3 = (ft, A,..., ft), E = (E , E ± , . . . , E d ), 
equipped with a symplectic structure 

d 

uj = ^2d(3i A dEi. 

i=0 

In the model example considered above d = 2, and we may put ft = 
1/T, ft = p/T, ft = -/i/T, and E = U,E 1 = V, E 2 = N. We iden- 
tify an abstract thermodynamic system with a Lagrangian manifold 
(possibly with a border): 

AcR 2(d+1) (ft£), 

and denote by i\ : A — > IR 2(ci+1 )(ft E) the canonical embedding of the 
system into the thermodynamic phase space (IR 2 ( d+1 )(ft E), oo). 
The conditions of A are as follows: 

(i) The manifold A is connected, simply connected, and admits 
E = (E , Ei, ... , E d ) as global coordinates. 

(ii) For any A > and any a G A, there exists a point b G A such 
that (p(b),E(b)) = (/3(a), XE(a)). 

Here we write ((3(a), E(a)) G R 2( - d+1 ^((3, E) for the coordinates of a 
point a G A acquired in the ambient space, and the notation XE is just 
(XE , XE ± , XE d ), for A G M, and E = (E , E x , . . . , E d ) G R d+1 (E). 
The point 6 G A corresponding to a given a G A and A > in the 
condition (ii) is unique, and we denote it b = Xa. The coordinates 
(3 G M d+1 (ft are termed the intensive coordinates of the system, and 
the coordinates E G M d+1 (i?) are termed its extensive coordinates. 

Proposition 1. Assume the conditions (i) and (ii) are satisfied. Then 
there exists a canonical choice of an action function S : A — >■ R ; 

d 

dS = i\(^^ dE ^ 

i=0 
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determined by the condition 

S(Xa) = XS(a), 

for any A > and a G A. 

Proof. Fix an arbitrary point oq G A. Any solution of the equation 
for S is of the form: S(a) := x + f "52 i=Q PidEi, where 7 C A is a 
path from oq to a, and x G R. Such a path always exists, since A is 
connected. The value of the integral does not change under continuous 
deformations of 7 since A is a Lagrangian manifold, and, furthermore, it 
does not depend on a choice of 7 since A is simply connected. Therefore 
the space of solutions of the equation for S is parametrized by x G R. 

Since A admits global coordinates E = (E , Ei, . . . , E d ), we can de- 
fine a function S = S(E) such that S(a) = S(E(a)), for any a G A. 
The condition (ii) implies that A(a) = /3j(Aa), i — 0, 1, . . . , d, for any 
a G A, and A > 0. Since ft(cs) = (dS(E)/dE i )\ E=E{a)l i = 0, 1, . . . , d, 
we obtain a system of equations 

^-(X- 1 S(XE)-S(E)) = 0, 

where i = 0, 1, . . . , d. It follows that 

S(XE) = XS(E) + c(X), 

where c(A) is a differentiable function, A > 0. If we take another 
real parameter [i > 0, then, applying the scaling twice, we obtain: 
on one hand, S(XfiE) = Xfj,S(E) + c(A/z), and, on the other hand, 
S(XfiE) = X[fiS(E) + c(/i)] + c(A). Therefore, we have an equation: 

c(Xfi) = Xc(fi) + c(X), (1) 

where A,/i > 0. Differentiating it with respect to /i and cancelling out 
the factor A > yields: c'(A/i) = c'(/z), i.e. d is a constant function. 
Then c(/i) = pp, + q, where p, q G R, and substituting it into the 
equation (TjQ) yields: pXfi + q = X[pfi + q] + (pX + q). It follows that 
X(jp + q) = 0, i.e. p = —q, and we have: S(XE) = XS(E) + q(—X + 1). 
Hence: 

S(XE)-q = X(S(E)-q). 

The function So(a) := S(E(a)) — q satisfies So(Xa) = ASo(a), where 
a G A, A > 0, and it is a unique function with such a property, since 
all other variants S(a) differ from So(a) by a constant. □ 

Definition 1. The real function S : A — > R satisfying the conditions of 
the proposition 1 is termed an entropy of an abstract thermodynamic 
system A. 
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Remark. We do not have a "shift of a reference point" relative to 
which we perceive the entropy of an abstract thermodynamic system. 
The entropy is canonically defined as long as we accept the extensivity 
of our theory (the condition (ii)). 

It is now the right moment to formulate the third law of thermody- 
namics, which is basically a condition on A at infinity: 

(iii) For any sequence of points {a n }^L C A, such that the inverse 
absolute temperature j3o(a n ) — > +00, while Ei(a n ) = Ei(ao), i = 
1, 2, . . . , d, as n — > 00, there exists a finite limit lim^oo S(a n ), 
and its value is the same for every such sequence. 

Remark. The fact that the limit in (iii) is actually the same has 
quite serious implications (the Nernst theorem). For any a G A, 
if 7 = {a(t)}ig[o,+oo[ C A is a continuous curve (a process) start- 
ing at a(0) = a and approaching the absolute zero of temperature, 
Po(a(t))~ 1 — > 0, as t — > +00, then it cannot be realized as a finite 
union 7 = U™T 1 7i, n G Z >0 , of isothermic and adiabatic parts (recall 
that 73 is termed isothermic iff /^(a)" 1 is constant along a G 7«, and 
7j is termed adiabatic iff S(a) is constant along a G 7$). Intuitively, 
and adiabatic process is a "swift" change of the values of extensive 
quantities, and the isothermic process corresponds to another extreme: 
it is very "slow". An informal statement of the Nernst theorem has a 
deep flavour of ancient Greek philosophy: it is impossible to reach the 
absolute zero in a finite number of steps. 

Let S = S(a) be the entropy on A, and let S = S(E) be a func- 
tion determined by S(E(a)) = S(a), for every a G A, (entropy as 
a function of extensive coordinates). Let / C {0,1, ... ,d}. Denote 
I := {0, 1, . . . , d}\I, and write: 

I = {i < % x < ■ ■ ■ < i n -!}, I = {i n < i n+1 < ■ ■■ < i d }, 

where n = |/| is the cardinality of /. For a phase space point (0, E), put 

Pi ■= (A » Pil, ■ ■ ■ , Pin-l) aIld E I '■= ( E inl E in+ll • • • » E i d) , aIld defille & 

projection 

tt 7 : R 2(d+1) R d+1 , {p, E) 1 — y (p T , Ej). 

Note that due to the condition (i) we have: dim^A) — d + 1, but the 
condition (ii) implies: dim^A) ^ d. 

Take a point ao G A and a subset / C {0,1, ... ,d}. Assume that 
there exists a neighbourhood U C A of the point ao G U , such that 
dim(7r/f/) = d+1. Assume that if a G U then \a G U , for every A > 0. 
Then we can consider (/3j, Ej) as local coordinates around a G A and 
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define a function S(j : tti(U) — > R as follows: 

S£(fr(a), = S(a) - ]T A(a)^(a), 

for every a G £7. Observe that Sjj((3i (a), A.E/(a)) = S[j((3i(\a), Ej(\a)), 
where A > 0, and therefore, the extensivity S(\a) = \S(a) implies: 

S I u (0 I ,\E I ) = \S[ r (0 I ,E I ), 

for any ((3i,Ej) e tti(U) and A > 0. The neighbourhood (7 C A is 
described by the equations: 

^ = -dSlfa, Ei)/dPi, Pj = dS^j, Ej)/dEj, 

where % e I, j e I, and (/?/, Ej) e ni(U). 
Consider a particular case / = [d], where 

[d] := {0,l,...,d-l}, 

and assume that Ed(a) > for all a e £/ (for example, .E^ is the volume 
of the system V). Then we obtain S$(/3[d], Ed) = EdS$(f3[d], 1), and 

The last equation can be perceived as a description of the projection 
7Tq(U) C M. d+1 ((3), dim7T0(f/) = d. The way to recover U from the 

projection 7r^(U) is given by the formulae: Ei/E d = —dSff(/3[d], l)/d/3i, 
i E [d]. 

3. Reduction of degrees of freedom 

If we have a collection of abstract thermodynamic systems A a C 
R 2(d a +i)^(a) ; ^ rfa e Z>Q) a = o, l,...,n- 1, then it is possi- 
ble to construct other thermodynamic systems out of it. The most 
simple case is where the systems do not interact with each other. 
Then together they form a system (the direct product) in a phase 
space M 2 ( d+1 )(/3, E), d = d + d± + • • • + <i„_i, with intensive coor- 
dinates (3 := ((3( \ . . . , /^ n-1 )), and extensive coordinates E := 
(E(°\ E^\ . . . , E^ 1 ^), described by a Lagrangian manifold 

A := A(°) x A« x ••• x A^ 1 ). 

The symplectic form is u := J2a=o J2t=o dPi*^ A dE\ a \ and the entropy 
S = S(a) on A is just a sum of the entropy functions: 

n-l 

s ( a ) = ^2S a (a a ), 

o=0 
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where a = (ao, a±, . . . , a n _i) G A, and S a is the entropy function on A Q , 
a — 0, 1, . . . , n — 1. 

If the subsystems A Q , ct = 0,l,...,n — 1, begin to interact, then 
we can perceive it as if some of the extensive degrees of freedom are 
being "released". Consider a linear non-degenerate transformation of 
extensive coordinates: 

d 
3=0 

where % = 0,1, ... ,d, and C = \\Cij\\fj =0 G GLa+i(M.). This lifts to a 

canonical transformation to = Ylt=o dfii A dEi = X^=o dfil A dE[, if we 
put 

d 
3=0 

where % — 0, 1, . . . , d. As long as the subsystems do not interact, we 
can hold the values of the extensive coordinates E' , E[, . . . , E' d . An 
interaction with respect to the last d — d' coordinates, ^ d' < d, 
means that E' d , +l , E' d , +2 , . . . , E' d are allowed to vary, so that the total 
system A ends up in a state with a maximum entropy (the second law 
of thermodynamics) . 

Write the entropy function S, a G A, in terms of extensive coor- 
dinates E'\ S(a) = S'(E'(a)), a G A. If we start in a point a G A, 
and release the coordinates E' d , +1 , E' d , +2 , . . . , E' d , then the final state 
b = b(a) after a relaxation to an equilibrium satisfies: 

E'(b(a)) = E[{a), (^fP) = °> 

V dE'j J E'=E>(b{a)) 

where % — 0, 1, . . . , d', and j — d' + 1, d' + 2, . . . , d. 

Since P-(a) = (dS' (E') / dE-)\ E f =E f^ , i = 0, 1, . . . , d, we observe that 
Pj(b(a)) = 0, j = d! + 1, d! + 2, . . . , d, and that the corresponding 
values of E'-{b{a)), j — d! + 1, d' + 2, . . . , d, are determined by E[{b{a)\ 
i = 0, 1, . . . , d' . This allows to define a "reduced" entropy S' = S'(E'), 
E' = (E' , E[,..., E' d ,) as follows: 

S'(E'(b(a))) = S'(E'(b(a))), 

for any a G A. The values intensive coordinates /^(6(a)), % — 0, 1, . . . , d! , 
can be expressed not in terms of function S' = S'(E'), but it terms of 
the reduced entropy S' = S'(E'): ^{b{a)) = (95'(^')/ 9£! j)U'=s'(6(o))- 
It follows that we obtain a Lagrangian manifold A' := R 2 ^' +1 ^(/3', E'), 
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f3' := (/3q, (3[, . . . , f3' d ,), described by the equations: 

= dS\E>)/dE>, 

where % = 0,1, ... ,d', and the symplectic structure is given by Co' : = 
Y2i=o A dE[. This manifold is nothing else but an image of 

M := {a G A | &(a) = 0, i = 0, 1, . . . , d'} 

under the canonical projection tt : R 2 ( d+1 ) (/?', E') ->■ R 2 ^ d ' +1 \p', E'). 

In terms of the original coordinates (/?, -E 1 ), the submanifold M C A 
can be perceived as follows. We have a family of planes 

L« := {£ G R d+1 1 £ + + • • ■ + 9? ) E A = p}, 

where p G R, and the coefficients g*- fc ^ := (q^\ q^\ . . . , q^) form a 
(d + l)-dimensional real unit vector, k = 0,1, ... ,d'. We assume that 
the intersections 

are (d — d')- dimensional, where (po,Pi, ■ ■ ■ ,Pd') € R d ' +1 . Put 
■^Po,pi,---,p d i • — {a G A | E(a) G Lp Q Pl t p dl } . 

For every a G A, we can compute E'(a), and then take the unique 
^ P0 (a),pi(a),...,p d ,(a) containing £(a). 

The planes L p fe \ p G R, = 0,1,..., d', define basically how the 
subsystems of A interact with each other. Once we "switch on" the 
interaction, the initial point a G A starts to move remaining in the 
submanifold £p ( a ),pi(a),...,p d /(a), until it finds a point b = b(a) of maximal 
entropy. The collection of points 6(a), where a varies over A (the image 
of the retraction a h-» 6(a)), yields M. 

Example 1. Consider two systems A^ and A^. with extensive coordi- 
nates (U( a \ V( a \ N^), a — 0,1, (internal energy, volume, number of 
moles). If the two systems are put in a thermal contact, then their in- 
ternal energies can change as follows: E^ — > E^ +x, E^ — > E^ —x, 
where x G R is a parameter. The other extensive parameters (volumes 
and numbers of moles) in this case remain fixed. 

Assume that there is another channel for a change of state: E^ — > 
E^+y, E^ — > E^ —y,N l — > — qy, where y G R is a parameter, 
and q > is a constant. Then the total picture is as follows: E^ — > 
+x + y, ->■ -(x + y), and ->■ JV« - qy, while the 
other extensive quantities N^°\ V^°\ and are fixed. Observe that 
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£"(0) -\- £(!) is a i so fixed. Define a linear transformation of extensive 
coordinates: 

(E' , E[, . . . , E' 5 ) := (E<® + EW,NW, V^,V^, (E<® - £ (1) )/2, 

The corresponding transformation of ((3^ a \p^ a \ Jl^), a — 0, 1 (inverse 
absolute temperature, pressure over temperature, and minus chemical 
potential over temperature) is of the shape: 

:= ((/^ + /3W)/2, j^p* 1 ', /? (0) - 

Reduce the thermodynamic system A = A^ ) x A^ with respect to the 
last two degrees of freedom in the transformed coordinates. A point 
a G A retracts to a point b = b(a) satisfying the conditions f3' 4 (b(a)) = 0, 
and f3' 5 (b(a)) = 0. We conclude: 

= 0, £ (1) = 0, 

in a state of thermodynamic equilibrium. This example mimics the 
black body radiation in a thermostat: the temperature of the radiation 
acquires the temperature of the thermostat, and its chemical potential 
vanishes. 

Example 2. Consider a system A consisting of three chemical sub- 
stances (A, B, and C) with an entropy function S = S(U, V, N , Aq, N 2 ), 
where U is the internal energy, V is volume, and Nq, Aq, A2 are the 
number of moles of the substances A, B, and C, respectively. Consider 
the following chemical reaction: 

A + B +± 2C. 

Note that the stoichiometric coefficients are well-defined already on 
the level of phenomenological thermodynamics. The number of moles 
can undergo the following changes: N — > N — x, Aq — >■ Aq — x, 
N2 — > N2 + 2x, where x G R is a parameter. We observe that the 
quantities D := Nq—Ni and A" := Ao+Aq+A^ are fixed. Introduce the 
extensive coordinates (E' Q , E[, . . . , E' A ) = (U,V, D, N, N 3 ) and reduce 
the system with respect to A3. We have: N = (N + D — N 3 )/2, and 
A"! = (N — D — N 3 )/2, so the intensive quantity /3 4 corresponding to 
E' A = A" 3 is of the shape: (3' 4 = — (pq +^2) /2 + At 3 , where /x , Jti, and /x 2 , 
are the intensive coordinates corresponding to Aq, N 2 , A3, respectively 
(the minus chemical potentials over the absolute temperature). If we 
start in a point a G A, then the final point b = b(a) satisfies /3' 4 (b(a)) = 
0, i.e. we obtain an equation 
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describing a 4-dimensional sub manifold M C A. In physical chemistry 
one obtains different characterizations of this surface depending an a 
chosen model. Consider, for instance, the following one: U(b(a)) = 
U(a), V(b(a)) = V(a), and the values of N (b(a)), Ni(b(a)), and 
N 2 (b(a)), are of the shape N (b(a)) = N (a) — x, Ni(b(a)) = Ni(a) — x, 
and N 2 (b(a)) = N 2 (a) + 2x, with x being a solution of the equation 

(iV 2 (a)+2:r) 2 
(iV (a) -x)(N 1 (a) - x) ' 1 J 

where if G M>o is the chemical equilibrium constant (the powers of the 
factors in the numerator and the denominator are the stoichiometric 
coefficients). In our case this is just a quadratic equation on x. If K — > 
+0, then in the limit we obtain: (2x + N 2 (a)) 2 = 0, i.e. x = —N 2 (a)/2, 
so the limit of the entropy (as a function S on the manifold A) in the 
final state b(a) is of the shape: 

S(b(a)) = S(U(a), V{a), N (a) + N 2 (a)/2, N^a) + N 2 (a)/2, 0), 

i.e. the chemical equilibrium is shifted completely to the left. On the 
other hand, if K — > +oo, then in the limit we have: (x — N (a))(x — 
Ni(a)) = 0. Since the quantities N (a) — x and A r i (a) — x must be 
non-negative, we obtain: x = min{A^o(a), A r i(a)}. It follows that 
N (b(a)) = 0, if N (a) < N^a), and that Ni(b(a)) = 0, if N {a) > 
Ni(a). The chemical equilibrium is shifted completely to the right. If 
N {a) = JVi(a) then 

S(b(a)) = S(U(a), V(a), 0, 0, N (a) + N^a) + N 2 (a)), 

i.e. A and B have reacted completely (since there was no excess of any 
of the substances) and everything has turned into C. 

Example 3. Consider a system A = A^ ^ x A^ which is a direct product 
of two copies of the system considered in the previous example. We 
have the extensive coordinates (U^ a \V( a \ N$ a \ N[ a \ N^) for each 
sybsystem a = 0,1. Assume first that the chemical reaction is to- 
tally suppressed and that the subsystems can exchange the internal 
energies and the chemical substances. Introduce the extensive coor- 
dinates U := U<® + t/ (1) , and N t := A^ (0) + , i = 0,1,2. The 
reduction with respect to the coordinates D := (U^ — U^)/2 and 
Di := (A^ 0) - iVf } )/ 2 , * = 0,1,2, yields a submanifold M C A de- 
scribed by the equations: (3^ — = 0, and Jlf^ — Jip — 0, % — 0, 1, 2, 
i.e. the temperatures and the chemical potentials of the substances be- 
come equal in the subsystems. If we take a G A, then after a relaxation 
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to equilibrium we obtain b(a) G M. The restriction of the entropy func- 
tion S\m written in terms of the coordinates (U, N , Ni, N 2 ) 
yields a generating function S = S(U, V (0) , V {1) , N , N h N 2 ) of the La- 
grangian manifold A of the reduced system. 

Switch on the chemical reaction A + B ^ 2C (the chemical equilib- 
rium constant is K). In accordance with formula ([2j) of the previous 
example, our point b(a) G M should retract now to a point c(a) G M 
on a sub manifold M C A described by the equation ^/(NqNi) = K. 

It is quite of interest to consider the following special case (the Gibbs 
paradox). Let a G A correspond to 

(4°\a),N[°\a),N 2 {0 \a)) = (n,0,0), 
(Nl ) 1 \a),Ni 1 \a),N 2 {1 \a)) = (0,n,0), 

where n G R >0 . Let V^(a) = V^(a) = v and U( \a) = U ( ^(a) = u, 
where u, v G M>o- Denote the entropy corresponding to the initial 
point a G A as S in (a). Denote the entropy corresponding to the final 
point c(a) G M as S out (a; K). We have: 

Sin{ a ) — S(u, v, n, 0, 0) + S(u, v, 0, n, 0), 

where S is the entropy as a function of extensive coordinates from the 
previous example. At the same time: 

lim S out (a; K) = S(2u, 2v , n, n, 0), 

and 

lim S out (a- K) = S(2u, 2v, 0, 0, 2n). 

K— >+oo 

The mixing entropy 

S mix (a; K) := ^(a; K) - S in (a) 

can vary depending on the parameter K. Suppose now than the chemi- 
cal substances A and B are very- very similar. For example, we measure 
this similarity in terms of molar masses M and Mi, respectively. Let 
the chemical reaction A + B ^ 2C consist in creating a substance C 
with a molar mass M 2 = (M +M 1 )/2. Put 

e := | M - M 1 \/M 2 . 

Let .fT = i^(e) smoothly depend on this parameter in such a way that 
K(e) <C 1, if s < So/2, and K{e) 3> 1, if e > £o, where e G lR>o- 

Fix M 2 = M G M>o, and denote the entropy function of the chemical 
substance with this molar mass as Sm = Sm(U, V, N), where U is the 
internal energy, V is the volume, and A^ is the number of moles. If 
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e < eo/2, then we do not "see" a difference between the substances A, 
B, and C, and therefore: 

S mix (a; K(e)) m S M (2u, 2v, 2n) - 2S M (u, v, n) = 0, 

due to the extensivity of the entropy function. On the other hand, 
if e > Eq, then we can clearly state that A, B, and C, are different, 
and there appears an observable jump S mix (a; K(e)) > 0. If we put 
S(2u, 2v, n, n, 0) = 2Sm(u, 2v, n), and take 



(the Suckur- Tetrode equation for N moles of a monoatomic gas in three 
dimensions), then in the region e > Eq we obtain: 

S mx (a;K(E))^2k B Nlog(2). 

This is precisely the jump of entropy considered in many discussions 
about the Gibbs paradox. The present example mimics a resolution 
of the Gibbs paradox in statistical physics suggested in [25]. This 
resolution is based on a construction of a kind of "number-theoretic" 
Bose gas of fractional dimension rf = 27 + l,0<7^1, but loosely 
speaking the philosophy can be reformulated as follows. Take a gas of 
identical "red" particles and a similar gas of identical "blue" particles 
(assume that all other parameters like mass, size, etc. of the particles 
are the same). Then the result of mixing is not a gas of particles 
some of which are "red" , and some of which are "blue" . The "correct" 
answer: the result is a gas of identical "purple" particles. These colours 
correspond to different dimensions of the number-theoretic Bose gas. 
In the present example, the substance A is, for instance, "red", the 
substance B is "blue" , and the substance C is "purple" . The chemical 
reaction A + B 2C can be perceived as a "loss of identity". 



5- 



4. QUASITHERMODYNAMIC FLUCTUATIONS 

Take an abstract thermodynamic system A C R 2 ( d+1 ) (/?, E) (a La- 
grangian manifold with respect to the canonical symplectic structure 
Sf=o dfiiAdEi satisfying the conditions (i), (ii), and (hi)). A reduction 
with respect to the extensive coordinates Ed'+i, £^'+2, • • • > Ed, where 
d' < d defines a retraction A — > M, a h-> b(a), to a submanifold M C A 
of dimension dim M = d' + 1. 

In quasithermodynamics the reduced thermodynamic quantities Ej, 
j = d! + 1, d' + 2, . . . , d, are not exactly fixed, but they fluctuate around 
the equilibrium values Ej(b), j = d' + 1, d' + 2, . . . , d, where a G M. 
The fluctuations are described by a collection of random variables 
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SE[d'+i,d] '■= (dEtf+i, 5Ed'+2, ■ ■ ■ , SEd). It is convenient to introduce 
the following concept. 

Definition 2. An abstract thermodynamic system A C M. 2 ( d+1>) (f3 , E) 
with an entropy function S = S(E) is termed linearly stable iff for 
every I C {0, 1, ... , d}, 1 ^ |/| ^ d, and oG A, the symmetric matrix 
obtained from the symmetric matrix \\(d 2 S(E)/dE i dEj)\ E=E ^\\ d :j=0 by 
deleting the rows i £ I and the columns j £ I , is negative. Denote this 
matrix as Sj(a). 

Assume that our thermodynamic system A is linearly stable. Then, 
in particular, S" d , +ld ^(a) < 0, a G A, where 

[d' + l,d] := {d' + l,d' + 2,...,d}. 

In quasithermodynamics the random vector 5E[ d >+i,d] is taken to be 
Gaussian with the joint density distribution function 

JSE [d , +1A {x,a) - (2irk B y d - d ')/ 2 * 

^ d-d' 

X 6XP { ~ 2k~ ^ x i[- S [d'+i,d\( a )]ij x j}> ( 3 ) 

where x = (xi,x 2 , ■ ■ ■ ,Xd-d>) £ R d_d ', and a G M C A. 

Remark. Note that by introducing a new physical constant k E into the 
theory (the Boltzmann constant) we actually make a step outside the 
paradigm of phenomenological thermodynamics. 

It is quite remarkable, that (according to Einstein) the associated 
intensive quantities @[d'+d\ '■= (Ai'+i, Pd'+2, ■ ■ ■ , fid) fluctuate as well. 
The equilibrium values in this case are (3\d'+i,d\{a) = (0, 0, . . . , 0), a G 
M C A, and the fluctuations Sj3[d'+i t d] '■= (SPd'+i, 8Pd'+2, ■ ■ ■ , 8(3d) form 
a Gaussian random vector with a joint density distribution function 

fw [d >+i, d] (y) := — !„ Md-dn/lt J„-i-f Oil /„MM/2 



(27r^)(^-^)/2(det[-^, +M (a)])V2 

d-d' 

GXP { " 2k~ ^ ^([-^'+i,d]( a )] _1 )«%}' ( 4 ) 



1,3=1 



where y = (y u y 2 , . . . , ya-d') d ' , and a G M C A. 

The classical probability theory (the axiomatics of A. N. Kolmogorov) 
identifies random variables with measurable functions on a space of 
events. More precisely, there is a probability model (f2, J 7 , P), where f2 
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is a set of elementary events, J 7 is a cr-algebra of events on ft, and P 
is a probability measure on (ft, J 7 ). The fluctuations and <5/3j, j = 
d' + 1, d' + 2, . . . , d, are measurable functions 5Ej : (ft, J 7 ) (R, B(R)) 
and 5^- : (ft, J 7 ) ->■ (R,S(R)), where 6(R) is the Borel a-algebra on 
the real line R. 

Why would we actually assume that the fluctuations 8Ej, 8{3j, j G 
[d 1 + 1, <i] should be modelled on a single probability space (ft, J 7 , P)? 
In quantum mechanics, if we consider the fluctuations of coordinates 
and momenta, this is not even possible (invoke Bell's inequalities 
and the Kochen-Specker type configurations (for an example, see [2"7])). 
There we have a statistical operator p = p> ^ 0, Trp = 1, acting on 
a Hilbert space H, in place of a probability measure P on (ft, J 7 ). In 
other words: 

It is natural to expect that if the thermodynamic system becomes 
smaller and smaller, then the probability model describing the 
fluctuations of extensive and intensive quantities becomes more 
and more quantum. 

More precisely, we expect that we can attach a Hilbert space H(a) 
to every point a G M C A, and that the fluctuations 5Ei and 5(3j 
can be represented by self-adjoint operators Q. t and Pj, respectively, 
i,j G [d' + l,d]. In the present paper I investigate the most natural 
possibility to define these operators. 

Theorem 1. Let A C M 2 ( d+1 ) be an abstract linearly stable thermo- 
dynamic system with an entropy function S = S(E), where E = 
(Eq, Ei, ... , Ed) are the extensive coordinates. Let M C A be the sub- 
manifold corresponding to a reduction of degrees of freedom associated 
with {Ej} j( z[ d / + i^, where d! < d. Put H(a) := L 2 (M. d ~ d '), for every 
a G M. Then there exists ijj(x; a) G T-L(a), where x = (xi, X2, ■ ■ ■ , Xd-d 1 ) 
varies over R d ~ d , such that 

fo K+M (^;fl) = \^(x;a)\ 2 , fsp [d , +1A {y;a) = \^(y,a)\ 2 , 

where ip{y;a), y = (yi,y 2 , ■ ■ ■ ,Vd-d>) € R d ~ d ' , is the (2k B )-Fourier 
transform of ip(x; a) with respect to x. 

Proof. Recall, that /i-Fourier transform of f(x) G L 2 (R n ) is defined as 
follows: 

(p{y) := {2nh) n/2 [ dx e- iyx/h (p(x), 



18 



A. E. RUUGE 



where h > 0, y G R n , n G Z >0 . If A = ||Aj||"j=i > is a constant 
matrix, then for 

, 2"/ 4 (detA) 1 /4 r 1 " . 

where x = (xi, X2, ■ ■ ■ , x n ) G R, we have: J* Rn c&c ^(x)! 2 = 1, and 
the /i- Fourier transform of y?h(:r; A) with respect to x is of the shape 
<fh(y, A) = (fh(y; Ar 1 ). Take a G M C A and substitute: 

n = d-d', h = 2k B , A = - S" d , +lyd] (a) . 

It is straightforward to check that 

V>(x; a) := y?2fc B (:r; --^, +M (a)) 

realizes the statement of the theorem. □ 

Remark. If we perceive the function constructed in the theorem in 
analogy with a coherent state in quantum mechanics, then we make a 
deep philosophical "discovery" : the Planck constant h corresponds to 
the Boltzmann constant ks multiplied by two. 

Do other thermodynamic wavefunctions tp(x;a) G L 2 (R.fr d '), a G 
M C A, actually make sense or is it just a fancy property of Gaussian 
exponents? At least it is natural to expect that the functions "similar" 
to (p2k B (x; — S'/ d , +1 d ](a)) do make sense and can describe a deviation 
of the state of the system from the thermodynamic equilibrium. Take 
x° G M. d ~ d ' and y° G M. d ~ d ' and consider a complex thermodynamic 
"wavefunction" 

ij x o, y o(x; a) := e^^tp^x - x°; -Sf, +M (a)). 

Note that the corresponding Weyl-Wigner function (2k b in place of h) 
is a Gaussian exponent concentrated in a point (x°, y°) G M. 2 ( d ~ d \x, y). 
Can thermodynamic "wavefunctions" be complex? Put 

- d 
Qj := x, Pj : = -z(2£; B ) — , 

for j = 1,2, . . . ,d — d' (self-adjoint operators on L 2 (R d_c! '(x)) corre- 
sponding to multiplication and derivation). Then for the moments of 
the fluctuations SEd'+j and 5(3d>+j in a state of thermodynamic equi- 
librium a G M C A we have: 

((5E d , +j ) m )= [ dxr o ,o(x;a),Qf^ (x;a), 

jRd-d' 

((Wd'+jD = [ dx^oM «)> ^o,o(x; a), 

JRd-d' 
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where m G Z>o, j = 1, 2, . . . , d — d', and the star denotes the com- 
plex conjugation. In particular: (SE d > + j) = and (5f3(i>+j) = 0, 
j = 1, 2, . . . , d — d'. If we replace ipo )0 (x;a) with ip x o ty o(x; a), then 
we obtain: 

(SE dl+j ) = x° p (Sfa+j) = yl 

where j = 1, 2, . . . , d — d' . 

Consider a model example. Let A = A^ x A* 1 ) be a thermody- 
namic system consisting of two subsystems with entropy functions 
S {a) = 5'(«)([/(«) ) y(«) ) AT(«)) ) a = 0,1 (the arguments are internal en- 
ergy, volume, and number of moles). Put (E , Ei, . . . , E 5 ) = (U^ + 
[/(i) ? v (°) , JV(°) , , , (U® -U^) /2) , and reduce with respect the 
coordinate E$. We have d' = 4. The systems are in thermal contact. 
Consider a discrete analogue of the process of the exchange of energy. 
Suppose that it takes place in "quanta" and in discrete steps in time: 
a subsystem a = can release a fixed amount of energy u > and the 
system a = 1 then absorbs it, or vice versa: the system a = 1 releases 
u > and the system a = absorbs it. If there is a tendency that one 
of the cases takes place more often than the other then we obtain a non- 
zero value of (SEs), which can be interpreted as an existence of a flow 
of internal energy from one system to another. The intensive coordi- 
nate p B corresponding to E 5 = (f/ (0) - £/ (1) )/2 is /3 5 = /3 (0) - where 
f3 {a) is the inverse absolute temperature in the subsystem a = 0, 1. A 
tendency to observe more often that one of the quantities (3^°\ (3^\ is 
greater than the other, leads to a non-zero value of (5(3 5 ). It can be 
interpreted as a gradient in inverse temperature, i.e. as an existence of 
a thermodynamic force. In a state of thermodynamic equilibrium the 
thermodynamic forces and flows vanish. We conclude that a complex 
thermodynamic wavefunction ip x o . y o(x; a) could describe a state near a 
thermodynamic equilibrium a G M C A with non-zero thermodynamic 
forces and flows. 

Remark. For the fluctuations 5(3j, 8E h j,l G [d 1 + 1, d], in the state 
^o,o (-; a) holds: 

j dxipQ (x; a)((PjQi + QiPj)/2)ipQfl{x] a) = 0, 

where j, I G [d' + 1, d}. Therefore, the linear correlation coefficient 

Corr^fij^Ei) = 0, 

for j, I G [d' + l,d]. This fact is compatible with the point of view 
of B. Mandelbrot [28] on the fluctuations of intensive thermodynamic 
quantities (for a review see [29|). Loosely speaking, one interprets the 



20 



A. E. RUUGE 



symbols 8/3 j, j G [df + l,d], as fluctuations of estimators of parameters 
of a probability distribution related to 8E h I G [d' + l,d], based on 
statistical samples. In Landau-Lifshits [30] one finds something com- 
pletely different: "their" fluctuations, which we denote A a /3j and A a Ei, 
j,l G [d' + 1, d], a G M C A, are linearly linked via the equations of 
state: 

A a (3j = £ A a ^, (5) 

dEjdE, E=E(a) 
l=d'+l 3 

where S = S(E) is the entropy as a function of E = (E , E%, . . . , Ed), 
j G [d' + l,d], and 

Corr(A a /3j, A a E t ) — 8 jtl , 

where 8jj is the Kronecker delta, j, I G [d r + 1, d]. We can now perceive 
this formula in analogy with semiclassical quantum mechanics, h — > 0. 
Having a semiclassical wavefunction ipn{q) = exp(zv4(g)//i)9?^(g), q G 
M n , where A(q) is a rea/ smooth function (the classical action), and 
<Pn(q) is a complex smooth function, it is possible to write it as a super- 
position of coherent states concentrated in the points of a Lagrangian 
manifold L C Mj^, p = (pi,p 2 , • • • ,Pn), 9 = (ft., 92, • • • , 9n), described 
by the equations: pj = dA(q)/dqj, j = 1, 2, . . . , n. The Landau-Lifshits 
equations are just a linearised analogue of these equations of classical 
mechanics. A generic quasithermodynamic "wavefunction" is a su- 
perposition of coherent states concentrated in different points (x°,y°), 
where x° = (x?, xjj, . . . , x° d _ d ,), y° = (yj, y§, . . . , y%_J), and and 
y° are values of A a E d i + i and A a /3d>+j, respectively, j, I G [1,<2 — d'], 
linked by the linear equations above (jSJ). A generic quasithermody- 
namic "mixed state" is a convex linear combination of orthogonal pro- 
jectors corresponding to quasithermodynamic "wavef unctions" . 



5. Tropical Pauli problem 

Speaking about a thermodynamic system A such as a one component 
gas S = S(U, V, N) (entropy S as a function of internal energy U , 
volume V, and number of moles N), we use the units of measurement 
which are adapted to the level of classical physics: length is measured, 
for instance, in centimeters, but not in angstroms or in parsecs, time is 
measured is seconds, but not in femtoseconds or billions of years. The 
number of moles in our system is measured in the scale of 1, 2, 3, ... , 
but we are not considering billionth fractions of moles or billions of 
moles. A convenient unit of measurement of absolute temperature on 
a classical level is one Kelvin. 
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On the other hand the approximate value of the Boltzmann constant 
in the CGS system is 

k B = 1.38 x W- 16 ergK-\ 

This numeric value is quite small and therefore the Gaussian expo- 
nents considered in the previous section are actually quite sharp. The 
standard deviations of the fluctuating quantities are proportional to 
y/ks, so the effect of the quasithermodynamic fluctuations is numer- 
ically in the 10~ 8 scale. Therefore we may formally consider k B as a 
small parameter, just like h is considered small in semiclassical quan- 
tum mechanics. 

Let us restrict to the case where we have only one reduced degree 
of freedom. In the notation of the previous section: d! — d — 1. For 
a fixed point a G A, we may also adjust the units of measurement of 
quasithermodynamic fluctuations 5Ed by taking a linear transformation 
— c(a)5Ed, c(a) G R>o, in such a way that the corresponding 
probability density function is of the shape: 

£<.)(a;) = (2irk B )- 1 / 2 exp(~x 2 /(2k B )), 

where xel Note that the physical dimension of the quantity is 
then [£( a )] = [k B ] 1 / 2 . Put rf a ^ := c(a)~ 1 S(3 d . Then the corresponding 
probability density is of the shape: 

f v(a) (y) = (2 7 rA; i? )- 1 / 2 exp(-y 2 /(2A; B )), 

where y G R. Note that = [/^b] 1 ^ 2 as well. The solution of the 

Pauli problem corresponding to the pair of functions (/^w,/,w) is of 
the shape: 

Vh (x) = (nh)-^eM-^/m), 
where one needs to substitute h = 2k B . We have: 

ftw(x) = \^2k B {x)\ 2 , f v (a)(y) = \^2k B {y)\ 2 , 

where x,y G R, and <fh is the /i-Fourier transform of <fh- 

In quantum mechanics, if we take a semiclassical wavefunction = 
f(x)exp(iS(x)/h), x G R, where / and S are smooth functions and S 
is real, then the H — > asymptotics of its /i-Fourier transform is (under 
some natural conditions) also a fast oscillating exponent. Assume that 
the equation y = dS{x)/dx has a unique solution x = x(y), so that we 
can define the Legendre transform 

S(y) = (S(x) -yx)\ x= x( y ), 

in every point y G R. Then for the h- Fourier transform we have ipniy) = 
exp(iS(y)/h)(g(y) + 0(h)), where g(y) is smooth. Let x(0) = and 
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write 

dS{x)/dx 2 = u{x), dS{y)/dy 2 = w{y), 
There is a link between the derivatives of the functions u and w: 

(^L-(ot)"=kL (6) 

where n G Zj>o- 

In tropical (older name - idempotent) mathematics, an analogue of 
the Fourier transform is the Legendre transform (the latter can be writ- 
ten as an idempotent integral). The link ([6]) mentioned was obtained 
for the real functions u and w, but we can formally apply it on complex 
functions u(x) and w(y). 

Observe that the quasithermodynamic "wavef unction" (p 2 k B {x) is 
of the form ^2k B { x ) = ( n ^ks)^ 1 ^ exp(iS(x) /(2k B )), where S(x) = 
ix 2 /2 is purely imaginary. For the (2fc^)-Fourier transform we have: 
& kB (y) = (n2k B )- 1 /*ex P (iS(y)/(2k B )), S(y) = iy 2 /2. Having this in 
mind, it is natural to define a tropical analogue of the Pauli problem as 
follows: given a pair of smooth real functions U(x) and W(y) defined 
in some neighbourhoods of x = and y = 0, respectively, find a pair 
of complex functions u(x) and w(y), such that 

Im(-u(x)) = U(x), lm(w(y)) = W(y), 

for which the formula (Q linking the derivatives in x = and y = 
holds. 

Example 4- A pair of constant functions u(x) = i and w(y) — i is a 
solution of the tropical Pauli problem corresponding to U (x) = 1 and 
W{y) = l. ' ' 

The insight is that one should reconstruct the functions not from 
the real, but from the imaginary parts of the function and its "Le- 
gendre transform" . In the introduction we have mentioned a somewhat 
counter intuitive fact about the Pauli problem: if it has a solution then 
there can be in fact several solutions. Let us look at what happens in 
the tropical case. Let 

u(x) = i( 1 + 2J — f[«m + «A m ]J , w(y) = ill + 22 — [w m + ip m ]j, 

where the real coefficients u m and w m are given, and the real coefficients 
A m and p m need to be reconstructed, m = 1,2, ... . The condition on 
the derivatives yields: 

d \ n 

i' l [w n + ip n 



i l u{x) dxJ i 1 u(x) 



x=0 
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where n G Z >0 , so the coefficients p n , n — 1, 2, . . . , become immediately 
known once we know A m , m = 1,2, ... . The first four equalities n = 
1,2,3,4 yield: 

—Wi = lm{ui + 
w 2 = Re{(w 2 + i\ 2 ) - 3(ui + «Ai) 2 }, 

w 3 = Im{(w 3 + i\ 3 ) - 10(ui + i\i)(u 2 + i\ 2 ) + 15(ui + iAi) 3 }, 

and 

- u> 4 = Re{(-u 4 + iA 4 ) - 15(ui + i\i)(u 3 + i\ 3 ) - 10(u 2 + i\ 2 ) 2 + 

+ 105(«i + i\if{u 2 + i\ 2 ) - 105(ui + iAx) 4 } 

From the first equation we find Ai = — w±, and substituting it in the 
second equation we obtain a condition: 

w 2 = u 2 - 3Re{(wi - iwi) 2 }. (7) 

From the third equation, taking into account Ai = —w 1 , we express A 3 
as a linear function of A2: 

A3 = IOM1A2 + w 3 — 10wiu 2 — 151m{(-ui — iwi) 3 }. 

Substituting this into the fourth equation, we obtain a quadratic equa- 
tion on A2: 

\j + 6wiWiA 2 + q = 0, (8) 

where 

q ■= ^~{ w 3 - 10u>iw 2 - 15Im[(ui - i^i) 3 ]} + j^i w 4 + u 4 - 

- 15u lU3 - I0u 2 2 + 105(> 2 - wl)u 2 - 105Re[(«i - m>i) 4 ]}. 

Since \ 2 should be real, there is an important condition on the discrim- 
inant: 

D := (Qu lWl ) 2 - Aq ^ 0. (9) 

It is of interest to notice, that if we express p 2 via Ai and X 2 , then we 
obtain: X 2 + (— p 2 ) = — Qu\Wi, i.e. — p 2 is just the second root of this 
quadratic equation. Denote the two roots mentioned A^, % = 1,2. 
Look now at the equations corresponding to higher degrees n = 

5, 6, If f(x) is a smooth function such that f(x) 7^ 0, then one can 

prove by induction (n ^ 5) the following formula: 



(— -V— = (-1) 

\f(x)dx) fix) { ' 



f (n) /7 (n-l) 



r3 f"f (n - 2) _ rp, (fTf {n ~ 2) , QnUj /',■■■, / ( - 3) ) 
°n+2 f n+3 ll n f n+4 + p n+l 
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where on the right-hand side we omit the argument x in the function / 
and its derivatives, Q n is a polynomial in n — 2 variables, C 2 +2 = (n + 
l)(n + 2)/2 and C^ +2 = n(n + l)(n + 2)/6 are binomial coefficients, and 
Tt n := C3, 2 = n{n+ l)(n + 2)(n + 3)/8 are the so-called tritriangular 

71 + 2 

numbers, the generating function: 

^ ^ m in 

/ J b rn X 



v ' m=l 



We have: 



(-D-^^ReK-ji 



(9 \2fc+l i 



9 \ 2fc +! 1 



x=0 



x=o' 



L -u(x) 9x7 i -1 ti(a;) 
so if fc ^ 2, we obtain from the first equation: 

(-l) fc ^2fe+l = Im{-[M 2fe+1 +iA 2fe+ i]+C' 2 2 fc+ 3[Mi+iAi](M2fc+^A2fc)+^2fe+l}, 

where Z 2 fc+i = Z 2k+1 (X 1 , A 2 , . . . , A 2 fc-i). This yields a linear link be- 
tween A 2 fc+i and \ 2k : 

A 2fc+ i - C 2fe+3 MiA 2fe + {(-l) fc w 2fc+ i - C 2fe+3 AiM 2fc + Im(Z 2fe+1 )} = 0. 

On the other hand, the second equation yields: 

= Re{-[w 2fe+2 + i\ 2k+2 ] + Cl k+4 [u x + i\i]{u 2k+ i + i\ 2k +i) + 
+ [ C 2k+i( u 2 + *A 2 ) - Tt 2k+2 (ui + iAi) 2 ](w 2fe + iA 2fc ) + Z 2fc+2 }, 

where Z 2fc+2 = Z 2k+2 (X\, A 2 , . . . , A 2 £_i). From here, invoking the linear 
link between A 2 fe + i and A 2 fc, we obtain: 

A 2 fe{Ai-u 1 [C 2A , +4 C 2fe+ 3 — 2Ti 2fe+2 ] + C| fc+4 A 2 } + Z 2k+1 = 0, 

where Z 2k+1 = Z 2k+ i(X 1: A 2 , . . . , X 2 k-i)- If the expression in the curly 
brackets in the last equation is not zero, then all the coefficients A 2 fc, 
and therefore A 2 fc+i, k ^ 2, are determined. It remains to notice that: 

C 2 2k+4 C 2 2k+3 - 2Tt 2k+2 = -2(k + l)(k + 2)(2k + 3), 

C| fc+4 = 2(fc + l)(A; + 2)(2A; + 3)/3, 

so we even obtain a condition — 3Ai«i + A 2 ^ which does not depend 
on the parameter k ^ 2. Taking into account that Ai = —w\ and that 
A 2 satisfies the quadratic equation A 2 + QuiW\\ 2 + q = 0, we reduce the 
condition just to q ^ (3wiU>i) 2 . 
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Assume the conditions ([HD, ©, and q ^ (3uiWi) 2 , are satisfied. We 
have the power series 1 + Ylm=i\ u m + iX m ]x m /m\ and 1 + J2m=ii w m + 
ip m ]y m /ml. To avoid a question about the convergence of these power 
series, it is natural to consider a truncated tropical Pauli problem. Let 
n G Z^o- Let U,W & C n °(R) be a pair of real- valued functions having 
a continuous derivative of order no in a neighbourhood of zero. Then 
the aim is to construct a pair of polynomials u(x) and w(y) over C of 
degree at most n , such that 

lm(u(x)) = U(x) + O(x no+1 ), lm(w(y)) = W(y) + O(y no+1 ), 

and the formula ([H]) induced by the Legendre transform linking the 
derivatives of order n holds for n ^ Uq. We term no the degree of 
truncation. 

Proposition 2. Let U(x) and W(y) be a pair of smooth real functions 
such thatU(x) = ^ + Y^=i u mX m / m\ andW(y) = i + J2m=i w ™y m / m - 
If the coefficients {M m }^ =1 and {w m }^ =1 satisfy the conditions: w 2 = 
u 2 — 3Re{(tii — iwi) 2 }, D > 0, and q ^ (3uiWi) 2 , where q and D are 
the constant term and the discriminant of the quadratic equation ([8]), 
respectively, then for any degree of truncation n G Zj> 2 , the truncated 
tropical Pauli problem has exactly two solutions. In the limit D — > +0, 
the two solutions coincide. 

Proof. See the explanations above. □ 

One may try to "improve" the formulae ([3]), (jlj), for the densities 
of distributions fsE d {. x ) anc ^ fspdiv) as follows. For a G M C A put 
S >2 (x;a) := S(E [0id _ 1} (a), E d (a)+x)-{S(E(a))+xdS(E)/dE d \ E=E(a) }, 
and $ >2 (y; a) ■= £>{<*} (#[0^-1] (a), (3 d (a) + y) - {S{d}{E[ 0jd -i\(a), d (a)) + 
ydS {d }(E [0id ^i](a),(3 d )/d(3 d \f] d= i3 d (a)}, where S^^d-i], /3 d ) is the Le- 
gendre transform of the entropy function S = S(E) with respect to the 
last argument, £[ ,d-i] = (E , E u . . . , E d _i), S{ d y(E[ 0id _i](a), d (a)) = 
S(E(a))— f3 d (a)E d (a), a G A (we assume that this transform is defined). 
Replace fsE d (x) and f S p d (y) with 

fsE d {x) = Z(a)~ 1 exp(S^ 2 (x;a)/k B ), 
fs/3 d (y) = ^(a)" 1 exp(-$ >2 (2/;a)/A;ij), 

where Z(a) and Z(a) are the normalizing factors. Note that in the point 
(x,y) = (0,0) we have d$^ 2 (y; a)/dy 2 \ y=0 dS >2 (x; a)/dx 2 \ x=0 = -1, 
and that is why we need a minus sign in the second exponent since 
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dS^ 2 (x; a)/dx 2 < 0, for x 7^ 0. More generally: 



/ d \ n + 2 ( 



1 d \ n 



y =o V S" 2 (x; a) dx ) S" 2 (x;a) 



x=0 



where S" 2 (x; a) := d 2 S^ 2 (x; a)/dx 2 , and n = 0, 1, We may assume 

that the coordinates and Pa are chosen in such a way that 

d 2 S >2 (x; a)/dx 2 = U(x), d 2 ^ 2 (y; a)/dy 2 = W(y), 

where the functions U(x) and W(y) are like in the proposition 2. Tak- 
ing n = 1 and n = 2 we obtain: w\ = u\ and w 2 = —u 2 + Suf. On the 
other hand, in the proposition 2 we have a condition w 2 — u 2 + 3[«f — 
= 0, so in the end this yields: 

w 2 = u 2 , 2u 2 = Zu\. 

This is a very strong condition on S = S(E) and there is no special 
reason for it to hold. Therefore typically one can not span a quasither- 
modynamic 2ks — > "wavef unction" over fsE d (x) and fsp d {y), and 
such fluctuations deviate from a pure state. 

6. Subtropical Pauli problem 

One may consider other improvements in the formulae for fsE d (%) 
and fsp d (y) corresponding to the fluctuations of the reduced degree 
of freedom (/3 d , E d ) in a linearly stable thermodynamic system A C 
R 2( - d+1 \f3 , E) . We keep the notation M C A for the submanifold cor- 
responding to (3d = 0. It is perhaps more natural to look for a thermo- 
dynamic "wavef unction" f 2 k B (x) of the form: 

if h {x) = exp(-x 2 /(2h))u h (x), 

where h = 2ks, and Uh(x) is smooth and expands in asymptotic power 
series in h — > 0, x e R. Assuming that Uh(x) does not grow too fast at 
infinity, we have a formula: 



■h d 2 
exp I — - 



u h (z) = (2nh) 1/2 / dx exp(-(x - z) 2 / \2h))u h {x) , 
Jr 



2 0z 2 

where zGl. Then for the h- Fourier transform (fh{y) we obtain: 

<Ph(y) = exp(-y 2 /(2h))w h (x), 
where Wh(y) expands into an asymptotic power series in h — > 0, 

'h d 2 



w h (y) ~ exp (__) Ufc (z) 



-w 
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where y G K. Let = exp(A/ l (a;)) and Wh{y) = exp(Bh(y)), where 



oo 



m=0 m=0 

where A m , / m , _B m , and g> TO , are real- valued functions, m = 0, 1, If 

these functions are analytical, 



n=0 n=0 



oc 



t-> I \ x * B m n y ^ \ g m ,ny 



n! 14 — ' n! 

n=0 ra=0 



then we can write: 

exp( J B /l (zx)) ~ exp ^- —J exp(A fc (a;)), (10) 

where This formula determines a link between the coefficients 

{4,n + «/m,n}m,n=o and { B m, tl + «? m ,n}«,„=o, and it suffices to assume 
that x varies in a small neighbourhood of zero. 

In the Pauli problem one assumes the knowledge of \(ph{%)\ 2 and 
\iph{y)\ 2 - It follows, that one can mimic it as follows: given the col- 
lections of coefficients {An,n}mn=o an d {B m ,4mn=oi reconstruct the 
collections {f m ,n}m,n=o and {9m,n} m,n=o- More precisely, given a pair of 
real- valued functions A h (x) and B h (y) admitting an asymptotic power 
series expansion in h with analytic coefficients, construct a pair of real- 
valued functions f h (x) and g h (y) which also admit asymptotic power 
series expansions in h with analytic coefficients, such that the formula 
(TinD holds for A h (x) := A h (x) + if h {x), B h {y) := B h (y) + ig h {y). 

We refer to the problem above as a subtropical Pauli problem (the 
small parameter h = 2k B describes a "deformation" of the tropical 
case). It is natural to consider also a truncated subtropical Pauli prob- 
lem in order to avoid the questions about convergence of the corre- 
sponding power series. Fix parameters mo, no G Zj>o (the degrees of 
truncation) and require that 

«* s »<«» " exp (55?) exp Aw) IL ~ °"* mo+1) ' 

for every n = 0,1,..., no- Basically it means that we consider the 
subtropical Pauli problem on polynomials. Note that without loss of 
generality we can also assume that / OT) o = and g m< o = 0, for all 
m = 0, 1, . . . , m , since the "wavef unctions" <Ph(x) and <Ph(y), h = 2k B , 
are defined up to a constant phase factor. One can also truncate only 
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with respect to h. In this case we require that a solution (f h ,g h ) is 
given by the polynomials in h of degree n ^ n , 

Let us look at what happens with the coefficients f m>n and g m ,n, 
m,n = 0,1,..., in more detail. Commuting the exponents on the 
right-hand side and then taking a logarithm in the equation (|TT|) yields: 

m=l 

Expanding the logarithm into a Taylor power series and then regroup- 
ing the terms, yields: 

1 /h\nJ^(-l)P+ 1 
I \ \ _ x 



x , ^ n 

B h {ix)^A h {x) + Y,- l {£fll 



n=l "" p=l P 

x "y n - r\\( d lM i d \ 2m \\ (12) 

2 -— ' mi!m 2 ! . . . mJ. L V dx dxJ J 

mi,m2,...,m p =l, ^ a=l 

mi+m,2H hm p =n 

Introduce now a notation. Let be a power series in x with real 
coefficients. Put: 

$ even (x) := + $(-x))/2, $ 0(M (x) := ($(z) - $(-x))/2. 

Observe that if x is real, then & even (ix) is real, while $ dd(ix) is purely 
imaginary. Separating the real and imaginary parts, we obtain: 

B h even (ix) + igfaiix) ~ + He(Q^(x)), 

-ifl*u(is) + ^(<x) ~ + lm(Q[ f \x)), 

where Qjf\x) is the series on the right-hand side in (|12|) . B^iix) ~ 
A/j(x) + Qjf (x). Since the series Qjf\x) is expressed via f h {x), we 
can conclude that once f h (x) is defined, we immediately know g h (x). 
Projecting on the even and odd parts, we conclude: 

B^ix) ~ ^Ln(^) + NQfW)U. (13) 

The coefficients at /i° yield: 

( J B ) e „ en (zx) = (A ) e „ e „(x), -i(B ) odd (ix) = (fo) dd(x). 

The first equality is a condition on the input data, and the second 
equality determines the odd part of fo(x). 
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Look at the coefficients at h 1 . The first equation in ( TT3l) together 
with (f )odd(x) = -i(B ) odd (ix) yields: 

1 ( d(f ) even (x )\ 2 

2 V dx 

where 



c(x), 



1 

+ 2 



(14) 



c(z) := -(.Bi) e „ en (iar) + (Ax) even (x)+ 

A Q (x)y | <9 2 A (a;) i lr dH(£o) 0(M (?a;)) i 2 

<9x / <9x 2 J even 2 1 dx 

Note that c{x) must vanish in a; = since (d(fo) even (x)/dx)\ x= o = 
0. This can be perceived as a condition on (Bi) even (0) — (Ai) even (0). 
Furthermore, there is an important condition: c(x) ^ 0. Since without 
loss of generality f h (x)\ x= o = 0, the function fo(x) is now known (if 
c(x) then there are two branches). The coefficients at h l in the 
second equation in (fT3l) determine (fi)odd{x): 

1. r/dL4 +i/o](a;) 



-i(Bi) dd(ix) 



(fi)odd(x) + -Im 



<9x 



+ 



- ocZcZ 



Look now at the coefficients at h k+1 for k ^ 1. The second equation 
in f lT3|) determines (/fc+Oodd^) since the coefficient in [lm(Q^\x))] odd 
at /i fc+1 depends only on {/z(x)}z<fc- The analysis of the first equation 
in (fl"3l) is slightly more complicated. In the definition of Q 



(/), 



a; we 



have a sum over n = 1, 2, 



. Let us extract the term n 
dA h (x) . d\\ 



1 explicitly: 



^ dx) 



a; 



where 



0(h 2 

(-^fc+l) evenix 



dx 

Then we obtain: 



:r 



(* + !)! 



2fc! 



<9 2 l4 fc + i/ fc ](a;) 
dx 2 



fe! 0[A, + z/,](x)0[A fc _ I + i/ Jfe _,](z) 



z=o 



l\{k-l)\ 



dx 



dx 



+ 



+ (Re(Q<' ) (s))) ( 



Therefore we obtain a linear equation on d(fk) e ven{x)/dx: 

d{f )even{x) d(f k ) even (x) f d(f ) od d (x) 8{f k ) odd (x) , 



<9x 

where Z^_/ X (x 



dx 



dx 



dx 



0, (15) 



requires a knowledge only of {/z(^)}z^fc-i- The left- 
hand side of this equality must vanish in x = 0, what is basically 
a condition on (B k+1 ) even (0) - (A k+1 ) even (0) , k > 1. If c"(x)\ x=0 ^ 
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0, then (d 2 (fo) even(x) / dx 2 )\ x= o 7^ 0, and differentiating the left-hand 
side of (fT5|) many times, one can find recursively all the derivatives of 
(fk)even(x) in x = 0. Without loss of generality, (f k )even(%) \x=o = 0. 

Theorem 2. Assume that (Bo) even (ix) = (Aq) even (x) , and that the 
function c(x) defined in (j!4p is non-negative in a neighbourhood of x = 
0, and c"(0) ^ 0. Then for any degrees of truncation (m ,n Q ) e (Z >0 ) 2 
there exists a polynomial P(h) in h = 2ks of degree m with constant 
real coefficients, such that the truncated subtropical Pauli problem for 
(A h + hP(h),B h ) has exactly two solutions vanishing in zero (polyno- 
mials in h = 2k b and x of degrees at most mo and n Q , respectively). 

Proof. See the explanations above. The polynomial P(h) is necessary 
to adjust the values of {B k ) even (0) - (A k ) even (0), k = 1,2,.... □ 

Example 5. Fix a point a G M C A, and consider a thermodynamic 
"wavef unction" of the shape: 

cp h (x; a) := (nh)~ l/4 exp(-(l + iha)x 2 /(2h)) 

where h = 2ks, and a = cr(a) is a real parameter. We have: <fh( x ] a ) — 
(7r/i) -1 / 4 exp(— x 2 /(2h)) exp(Ah(x)), where Ah(x) = —iax 2 /2. Then 
<p(y;a) = (7r/i)~ 1//4 exp(— y 2 /{2h)) exp(B h (y)) (the /i-Fourier transform 
of (Ph(—'i a )), where 

B h (y) = — log(l + iha) + . 

Kyj 2 sv J 2 1 + iha 

This yields the following data for a truncated subtropical Pauli prob- 
lem: 

A (x) = 0, A 1 (x)=0, B (y)=0, B 1 (y) = a 2 y 2 /2. 

For the quantity c(x) we obtain: c(x) = o~ 2 x 2 /2, so df (x)/dx = ±ax. 
Assuming that /o(0) = 0, we obtain: fo(x) = ±ax 2 /2. Put 

^\x; a) := (nhy^exp ( - ^[1 + (-l) a ik(a)]) 

where a = 0, 1, Look at the 2-dimensional subspace L = L(a) in TL = 
L 2 (M.) spanned over (pf\—;a) and (f^\—;a). Let us leave out the 
fixed point a £ M C A in the notation and write just ^ a \x) in place 
of ip^(x;a). The inner products between the basis vectors ip^\x), 
a = 0, 1, are of the shape: 

(ri°U 0) ) = i, (<p£\<pP) = i, (^\A l) ) = (i-^ 2 , 
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where one should take the branch (1 — iha) — > 1, as ha — > 0, of the 
square root. Consider an orthonormal basis in L = L(a): 

(1)/ \ (0)/ w (0) 

= *,?>(*), >(*) := V " (X) -\ (X) ^^ \ 

where we leave out the argument a in as well as in ip^\ ■ 

A generic statistical operator p concentrated on L = L(a) is (in the 
Dirac notation) of the shape: 



ra,n=0,l 



where the coefficients p m>n = p m , n (a) G C form a 2 x 2 matrix P = 
|bm,n||m,n=o,i such that P = P^ ^ and Tr(P) = 1. These additional 
data provided by the matrix P = P(a) define a "deformation" of the 
quasithermodynamic fluctuation theory in the point a G M C A. The 
variances ((5Ed) 2 ) and {(5(3d) 2 ) are computed as 



({5E d f)= E jw^SLio 2 !^ 

m,n=0,l 

«w 2 >= E p^n^SLiP 2 !^ 



m,n=0,l 

where Q = x (multiplication by x), and P = —i{2kB)d/dx. The 
covariance coefficient is computed as 

(5f3 d 5E d )= e PrnMZmttl), 

m,n=0,l 

where R := (PQ + QP)/2. If the statistical operator p — \ip^ B ) I' 
then we recover in the leading term in 2/c# — > the formulae of the 
standard quasit her mo dynamics. 

Note that the case c(x) = is actually admissible, but it is a little 
technical and we do not go into its details. A generalization of the trop- 
ical and subtropical Pauli problems to many dimensions is straightfor- 
ward, but their analysis involves some overdetermined linear systems of 
equations. In principle, the general picture where an existing solution 
can have an "antipode" remains valid. 

7. Conclusions and discussion 

The main motive of the present paper is that the probability model 
describing quasithermodynamic fluctuations becomes more and more 
" quantum" if the thermodynamic system gets smaller and smaller. By 
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that it is not implied that the Planck constant H starts to play a more 
visible role, but rather that the very "nature" of the probability model 
becomes different: 

There is no particular reason to perceive the fluctuations of 
thermodynamic quantities as random variables sharing the same 
probability space. 

It is natural to expect an analogy with quantum mechanics, i.e. that 
at some point one needs a statistical operator on a Hilbert space to 
describe these fluctuations. A more detailed analysis implemented in 
the present paper shows that the role of h can be taken over by the 
Boltzmann constant ks multiplied by two. 

We look at an abstract thermodynamic system A and reduce it to A' 
with respect to several degrees of freedom n = dim(A) — dim(A'). This 
leaves us a "heritage" consisting of fluctuating intensive and extensive 
thermodynamic quantities associated with these degrees of freedom, 
for which we state an analogue of the Pauli problem: reconstruct a 
"wavefunction" (more generally - span a class of statistical operators) 
consistent with experimentally available data (the marginal probabil- 
ity densities of the collections of quantities which can be measured 
simultaneously). 

Many important features of the thermodynamic Pauli problem can 
already be seen if n — 1 and in the quasithermodynamic ("semiclas- 
sical") limit 2ks — > 0. We formalize these observations in a theorem 
about a truncated subtropical Pauli problem proven in the main text. 
An important conclusion here is the following. In principle, to describe 
the fluctuations of thermodynamic quantities in analogy with quan- 
tum mechanics we need to attach a Hilbert space H{a) ~ L 2 (IR n ) to 
every point a G A' and look at statistical operators on this space. It 
turns out that if we are only interested in the first corrections to the 
standard quasithermodynamic fluctuation theory, then one can signif- 
icantly simplify things considering instead of the infinite dimensional 
Hilbert space %{a) a 2-dimensional subspace L(a) C H(a). Typically 
(i.e. under the conditions of the theorem mentioned) the truncated 
subtropical Pauli problem is going to have two solutions which corre- 
spond to a pair of quasithermodynamic "wavefunctions" tp^ (x; a) and 

t P2k B ( x '> a )' x e ^" Knowing one of the functions allows to reconstruct 
the second one explicitlly. The space L(a) is then 

L(a) := spaji c {cp { ^ B {-;a),cp { ^ B (--a)}. 

Instead of a generic statistical operator on T-L(a) one may approximate 
the fluctuations by a statistical operator p(a) concentrated on L(a) 
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(the operator p(a) in the latter case is, loosely speaking, nothing more 
but a 2 x 2 matrix). In the leading degree in particular cases one 
recovers the standard quasithermodynamic fluctuation theory, but the 
general formulae contain new interference terms, which are of interest, 
for example, in the quantum information science. 
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